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Abstract: In this paper we reconstruct convection coefficients from boundary 
measurements. We reduce the Beals and Coifman formahsm from a hnear first 
order system to a formahsm for the 9-equation. 

1 Introduction 

The pioneering work of Nachman and Ablowitz jEI) Sylvester and Uhlmann [2^, Nachman |17jand 
Henkin and Novikov ^U] introduced inverse scattering methods to the parameter identification 
problems. In their work, the linear Schrodinger equation in the physical space is paired with a 
pseudo-analytic equation in the complex space of the parameter. Another method, due to Beals 
and Coifman 2 , pairs a first order d system in the physical space with a pseudo-analytic matrix 
equation in the parameter space. Sung analyzed lower regularity assumptions in |22| \2\\ \ 1^ . This 
method was ingeniously used by Brown and Uhlmann |4j in unique identification of the conductivity 
o" in V • cjVu = and by Cheng and Yamamoto 0! HI in proving unique determination of the 
convection coefficients hi and 62 in Au -|- hiUx + b2Uy = 0. 

We consider here the scattering problem for d- equations (theorems 1.1 and 1.2 below). Here 
d = {dx + idy)/2 is the Cauchy-Riemann operator. This can be seen as a diagonal version of 
the formalism in Beals and Coifman, see lemma l2Tl Due to the symmetry between the scattered 
solutions in the physical space and the ones in the parameter space, we are able to present a 
non- linear analog of the Fourier inversion formula (compare and below). 

As an application, we revisit the inverse problem proposed in |6^ and present a reconstruction 
procedure. The method is based on solving a singular boundary integral equation in the Hardy 
space of functions in the exterior of the disc. This method was first introduced by Knudsen and 
Tamasan in connection with the electrical impedance tomography problem in The method 

presented here can be seen as its generalization. 

I was informed recently about the reconstruction step being obtained independently by Tong, 
Cheng and Yamamoto 0. I thank them for letting me know about their new result. The main 
difference of the method presented here from their method is the formalism of in inverse scattering. 

For /c G C arbitrarily fixed, we say that u behaves like e^^^ (written u ^ e^^^) in L'^(M.'^) for large 
z, if ti(z, /c)e~*^^ — 1 € L^{M.l). We use the notation (k) = (1 + |fep)-^/^. The scattering method is 
the content of the following two theorems. 

*This work was done during the author's visit at IPAM-UCLA in the Fall 2003 
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Theorem 1.1 (Forward Scattering). Assume that q G Lc(M^), p > 2 has compact support. For 
each k E C, the equation 

^{z) + q{z)^{z)=0, zee, (1) 

has unique solutions '^r{z,k) ~ e^^^ and '^i{z,k) ~ ie^^'^ in LP for large z, and the scattering 
transform 

t{k) = -- I e^q{z) {^r{z, k) - i^i{z, k)) dfi{z) (2) 

is well defined. Moreover, if q e Wc'^\^l) for some e > and k eC — {0}, we have 

W^riz, k)e-''>' - + W^iiz, k)]e-'''' - i||ip(K2) < C{k)-' (3) 

and 

mr{z,k)-i^i{z,k)]e-''''-2\\w,,,(j^2^ < C{k)-', (4) 

and then t G ^^(M^) for each r>2/{e + l). In paHicular t G L''(M|) n ^'■'(M^) n L^(M^) for some 
r < 2, where r~^ = r~^ — 1/2 and r'~^ + r~^ = 1. 

Theorem 1.2 (Inverse Scattering). Letq, ^r, *i andt{k) andr, r' , f he as given in the forward 
scattering. Then the equation 

^(fc)+t(fe)¥(fc)=o, fceC, (5) 



has unique solutions $r ~ e*^*^ and $j ^ ie^^^ in L'^(M|) for large A; G C. Moreover, ^'s and 
are related by 

Re^>i = -Im^'r, Re^>r = Re^r, (6) 

m particular $r ~ = ~ a'^^^ 

= -- /" e'^i(fc) (^r(^, A:) - i^i{z, k)) dfi{k). (7) 

Let C be a bounded, simple connected domain with Lipschitz boundary and p > 2. For 
61,62 G LP{n) and c/ G W^-^/P'P{dn), let u G H^^.p be the unique solution of the boundary value 
problem 

An(ar) + 61 (ar) — (ar) + 62 (a;) ^ (x) = 0, x G O (8) 
0x1 0x2 

u{x) = g{x), X G dQ. 
The Dirichlet to Neumann map A^^^b^ : W'^-'^/P'P{dn) W'^-'^/P'P{dn) is given by 

^hub29{x) = i'i{x) — {x) + f2(x) — (x), X G dn, 
0x1 0x2 

where (1^1 (x), ^2(3^)) is the outer normal at x on the boundary. Cheng and Yamamoto proved that 
Afej uniquely determines 61 and 62 in LP{p,). 

Working with the equation in the whole plane and using the inverse scattering for (^-equations 
allows us to go beyond uniqueness and present a method of reconstruction. We prove the following 
result. 
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Theorem 1.3. Let C &e hounded, simple connected domain with Lipschitz boundary, and 
61,62 £ Wc'^i^), with support inside 0, for some e > . Then 61,62 can be reconstructed from 
^61,62- 

The fact that they vanish on the boundary is not a severe restriction as one can always extended 
the coefficients across the boundary, preserving the regularity, and then have them vanish outside a 
ball. The Dirichlet-to-Neumann map can be pushed to an outside boundary as shown by Nachman 
in see also ^S]- While L^(0) is enough regularity to prove unique determination of 61,62, we 
assume here e-extra regularity and provide a reconstruction method. 

In the end we point out the connection with the first order d system and characterize its Cauchy 
data in terms of the Dirichlet-to-Neumann map of a related second order elliptic equation, thus 
answering a question of Uhlmann in j25j . 



2 Proof of the theorems 11.11 and 11.21 
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We identify a point in with a point in the complex plane by xi + ix2 = z. By d " we denote 
the solid Cauchy transform 

rV(^) = - / ^d^(c), (9) 

where is the Lebesgue area. We also denote by e{z, k) = exp{i{zk + zk)). 

We look for solutions of (QJ of the form ^'^ = V're*^^ and = ii/^je*^^ with ipr, i^i € 1 + LP(M^). 
The equations for ipr respectively ipi are 

d 

—Tpr + qe(z, -k)'ipr = 0, (10) 

d — 

—Ipi - qe{z, -k)ipi = 0. 

A key ingredient is the Hardy-Littlewood Sobolev inequality which yields d ^ : L^(M^) 
is bounded (see Stein ^0 ) for p and p related by 

111 

-~ = --l^- (11) 
p p 2 

Since q G Lf(R2) c L2(]R2) and LP(R2) • L^{R'^) c LP(R2) we have d~^{q-) : LPt{R'^) LP{R'^) is 
a bounded operator. Since q has compact support we can use Rellich imbedding to conclude that 
d~ {q-) : LP{R^) L^iR^) is compact. Then we can apply Predholm's alternative in L^iR?) to the 
equivalent integral equation 

{I + d-'[qi-)e{; -k)-mM^) - 1) = d-'[qei; -k)], (12) 
{I-d-'[qi-)ei;-k)VmM^) - 1) = -d-'[qei;-k)]. 

The fact that the homogeneous equation has only the null solution comes from Liouville's theorem 
for pseudo-analytic functions with coefficients in L^(R^) nL^(M^) shown by Vekua |V62j . Since we 
integrate in Q over the support of q, together with the imbedding Lf^^ C Lf^^, gives a pointwise 
well defined t{k). 
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For A; € C — let {d — ik) ^ be defined by e(z, k)d (e(z, —k)-) and let the indexes p and p be 
related by An interpolation (with e being the interpolation parameter) between the estimates 
of Nachman m W - ik)-' fWiP < ^^II/IIlp(r2) and \\(d - ik)-^ fW^P < (C/|A:|)||/||h.i,p(m2) gives 

\\{d-ik)-^f\y<^J\f\\^,,,^^,y (13) 

See Proposition 2.3 in ^31 for details. Since |e(2;, A;)| = 1 the last estimate implies that 

\\d ^{e{-,-k)q)\\LP(R2) <C{k)~'\\q\\we,p(^R2y 

The decay rate in ^ follows from the (uniform in k) bounded-ness of the map [I — d ^qe{-, —k)-]~^ 
from L^(R^) to L^(M^) as explained above. The further regularity property for the combination 
+ in ® will be shown in Lemma l2. II For now we assumed it holds. 

Brown and Uhlmann ^ showed that q E Lc(IR^) implies t G L^(R|). While this is good enough 
for existence, for reconstruction we need t ^ U for some index r < 2. This is ensured by extra 
regularity imposed in q as was shown by Knudsen and the author in ^31- For completeness we 
repeat the arguments. The main ingredient is an LP' bounded-ness property for pseudo-differential 
operators with non smooth symbol (see Coifman and Meyer [H] or Brown and Uhlmann (2|)- If Mq 
is defined by 

Mq{k) = — / e{z,k)q{z)a{z,k)dfi{z), 



vr _ 

where a has compact support in z and ||a(-, A;)||j|^i(]g2) < Cik)^'', then M : L'^{M.1) L^(M^) is 
bounded. 

Rewrite now 

tik) = -2iTiq)i-2ki,2k2) + Tik) (14) 
T{k) = - [ eiz,k)qiz)[Mz,k)+iJi{z,k)-2]dfiiz), (15) 

where T is the Fourier transform. Since q G if (M^) C L''{RI) for 1 < r < 2 we get J^{q) G L''(M^) 
for s > 2/{l + e). 

Let M be the operator defined by a{z, k) = x{z)[i>r{z, k) + ipiiz, k)-2]e Wc '^(M^) C H^{R'^), 
where x is a cut-off function equal to 1 on the support of q. The following chain of inequalities for 
< 6 < e give the result 

II^IIl=(m2) = Mq\\ ^s(^r2) < C'||(A:)"^||^(i/,_i/2)-i ||g||i2(K2), 

for 6 > 2{l/s - 1/2) or equivalently s > 2/(5 + 1) > 2/(e + 1). 

In order to exhibit the relation with the old formalism, we prove theorem II. 21 bv reducing it to 
the former. Let us define mi{z, k) and m2{z, k) in terms of the ip^s by 

mi{z,k) = ^{ijr{z,k) +iJi{z,k)) (16) 
7712(2:, A;) = ^e(z, -k) {ipiiz, k) - ■^^(z, k)) 

The simple result below shows that (7711,7712)* is the first column of the Jost matrix in the complex 
geometrical optic solutions of Beals and Coifman. 
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Lemma 2.1. Let mi and m2 defined in ()16() . Then mi{-, k) — 1, m2{-, k) E LP(M^), and they satisfy 

dnii = qm2 (17) 
{d + ik)m2 = gmi. 

Moreover, the following estimates hold, 

\\mi{-,k)-l\\wi,P(^^l)<C{k)-' (18) 
\\m2{;k)\\LPiRl)<C{k)-\ (19) 

Proof. From their definition mi(-, A;) — l,m2(-, fc) € L^(IR^) since 

V'r,V'i e -^^^(IK^) and \e{z,k)\ = 1. 
The fact that they solve the system (|17|) comes from a straightforward calculation and the equations 
(|10)) . The estimates of decay in k for noth mi and come from the estimates ((3| for ^/^^ 
and ipi proven above. We are left to justify the extra smoothness gained by mi. From the first 
equation we have that mi — 1 = 5 {qm2). Since q G Lc(]R^) C L^(M^) and L? ■ C LP we have 
d {qm2) G ^^^'^(M^) with an imbedding constant which depends on the support of q but it is 
independent of k. We have the following chain of inequalities. 

\\mi{-,k) - l\\w^,P(R2) = \\d ^{qm2)\\w-^,P(M.2) < C\\q{-)m2{-,k)\\Lp(^2j 
< C|kllL2(iR2)||m2(-,A;)||ip(K2) <C <k>-' . 

□ 

This also completes the proof the theorem ll.il 

Formulate the inverse scattering formalism of Beals and Coifman only in terms of the first 
column of Jost matrix, see Knudsen and Tamasan for details. For the analysis with q € 
L1(]R2) n L°°(M2) see Sung or Brown and Uhlmann [1| for g G (M^). 

Theorem 2.2 (Beals & Coifman scattering method). Let q G Wc'^{M?), For any z ^ C, 
the system ()17() has a unique solution m.i{z,k), m2{z,k) with {mi{-,k) — l,m2{-, k)) G L^{M.1). 
Furthermore, the map k m{-,k) is differentiahle (in the norm topology) with values in W'''^{dx) 
and satisfies pointwise in z £ C the system 

d 

^mi(2;, k) = t{k)e{z, -k)m2{z, k), (20) 

d 

-^m2{z, k) = t{k)e{z, —k)mi{z, k), 
dk 

where 

t{k) = -- j e{z,k)q{z)mi{z,k)dn{z). (21) 

Jk.2 

Look for solutions of © in the form ^{z,k) = ie'^^^(f)r{z,k) respectively = e'^^^(f)r{z,k). As 
in the forward problem, they must satisfy an integral formulation analogous to (|12|) where the role 
of k and z is reversed. Since t{k) G L'"(]R|) fl L'^(M|) we have existence and uniqueness for their 
solution in L''(]R|), where = r^^ — 1/2. Using the equations (|2flj) it is easy to check that 

d 

^(mi — m2){z, k) = —t{k)e{z, —k)mi — m2{z, k) (22) 
dk 

d 

-T^imi + m2){z, k) = t{k)e{z, -k)mi + m2{z, k). 
ok 
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By the uniqueness result we must have 



(j)i(z,k) = mi(z,k) + m2{z,k) (23) 
(priz, k) = mi{z, k) — m2{z, k). 

The following equalities show the relation between solutions of the forward and inverse equation. 

$. = ie"^(t)i = ie*^'=(mi + 7712) = ^"(^r + ^i) + -^^i^^ -k){Wi - Vv) 
i 1 1 i 

Similarly, = Re^'r — iRe^'j. These prove the identities Formula ^ is due to a symmetry 
argument as follows. Starting with q produce ipr and ipi by solving (0). Via © produce 0r and 
and then as in ((T)). Take this t(A;) and do now forward scattering starting from the A;-space, i.e. 
produce and by solving © and via (jH)) produce "ifr and ^'j. Define a potential gi using 
for the z-space. In particular we know that for any A; G C we have d^r + Q^r = since we started 
that way, but also now we have d"^r + Qi^r = 0. In particular we have {q{z) — qi{z))"^r{z, k) = 
for all /c € C. Hence q = qi- 



3 Reconstructing convection coefficients 

In this section we apply the above scattering method to reconstruction of the convection coefficients 
h , 62 in 

Au{x) + bi — {x) + b2{x) — {x) = 0,xen (24) 
ox ox 

from the Dirichlet-to-Neumann map A^^ fe^- Here C is a bounded, simply connected domain 
with Lipschitz boundary. 

We assume here that 61, &2 £ Wc'^{^), p > 2 are real valued maps with compact support in 0, 
and set b = {bi + i62)/4. 

The following result from Vekua |V62j makes the reduction of (|24j) to a 9-equation. If u is a 
solution of (|24() then w = du solves 

dw{z) + b{z)w{z) + b{z)w{z) = 0. (25) 

Lemma 3.1. Let Q be simply connected with Lipschitz boundary. If u G VF^'^(r2) is a solution of 
(|24j) ■ then w = du £ W^'^{i^) is a solution of (|25|1 . Conversely, if w £ W^'P{^}) is a solution of 
H25|) then there exists an u £ W'^'^{Q) solution of (|24() and such that du = w in J7. 

Proof. By Sobolev imbedding we have u £ C^+"(r2) with a = 1 — 1/p and w £ C"(J7). As a 
direct consequence of the Poincare lemma, notice that if dw is real valued, then w = du for some 
real valued u. Indeed 2dw = {dx + idy){f + ig) = {dxf — dyg) + i{dxg + dyf). By assumption 
dxg = —dyf, from where the one- form gdy — fdx is exact. Therefore, there exists a real valued 
F such that dF = {-f)dx + gdy. We have w = f + ig = dx{-u) - idy{-u) = d{-2F). The 
equivalence is now apparent. □ 

Now we extend b £ Wc'^{0,) by zero outside O. Its extension denoted also by b preserves 
regularity b £ Wc'^{M.'^). From now on we shall work with solutions of H25|) in the whole plane. 
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Lemma 3.2. The equation ()25() has unique solutions in the whole plane Wr{z,k) ~ e^^^ respec- 
tively Wi{z,k) ~ ie*^*^ in LP{RI) for large z. Moreover, e'^^'^Wr — l,e~*^^Wj — i G W^'P{R1) and 

Wr{;k),W,{;k)eWl;P{R''). 

Proof. As in the proof of theorem 1.1, we look for solutions W{z, k) = e^^^w{z, k) with w — 1 G . 
The equation for w is 

d{w{z) - 1) + b{z){w{z) - 1) + e(z, -k)h{z){w{z) - 1) = -b{z) - e(z, -/c)6. (26) 

Using the fact that d ^ : f £ (M^) i-^ W^''p{M?) together with b of compact support we get 
d (b-) : L^(R^) L^(M^) is a compact operator. We apply Fredholm's alternative in L^(M^) to 
the equivalent integral equation 

{[I + d-\b{-) + e{; -k)m}iH^) - 1) = -d~\b + e{; -k)b]. 

Uniqueness comes from Liouville's theorem for the (^-equation with coefficients in L^(M^) n-L^(M^), 
see Vekua |V62j . By construction we already have that g = Wr — 1 £ W^'P{W1). Then Wr{z,k) = 
e''^{g+l) G dWr = ike''^g{z, k)+ike^''' +e^'^dg G Lf^^CM^) and = e*^% G L^i^^i^D- 

Similar relations hold for Wi. □ 

To simplify notations, let 

q{z) = 6(z)e^"'''(")-^"'^(") (27) 

denote a new potential and notice that if t(; is a solution of H25|) then v = e^ is a, solution of 

dv + qv = 0. (28) 

Since b G LP[m?) n LP(M2) we have that d~^b G L°°(M2) n C^, see Vekua |V62] . Then e"^'"'' G 
L°°(M2) and so g G LP(M2) n U>{R^) . 

The next theorem relates scattering solutions of (|25|) to scattering solutions of (|28j) and gives 
the behavior in /c of Wr{z, k) and Wi{z, k). 

Proposition 3.3. Let b G Tyc'^(IK^); /or some e > 0. Let W,. and Wi be the scattering solutions for 
(|25|) as given by the lemma above, and let o-nd he the scattering solutions of (|28|) as given 



by the theorem 1.1. Then Wr = e~^ ^"^r, Wi = e~~^ ^^i and 

\\Wr{z,k)e-^^^ - e-^"li,(K2) + \\W,{z,k)e-^^>^ - ie~^~\LPm) ^ C'(A:)^^ (29) 
[Wriz, k) - iWi(z, k)]e-'''' - 2e-^''^ < CikY 



Proof. The fact that Wr and Wi solve (|25p is trivial. Uniqueness result of lemma EIH ensures that 
they are the scattering solutions of ()25() . The estimates follow directly from the estimates for ipr 

and ipi in and Q and from the fact that ^ G L°°{M.1) as noticed before. Again, the 

imbedding W^'P{R^) C C^-^/p(M?) shows that the estimates ^ hold pointwise in z G C. □ 

We have now all the ingredients necessary for reconstruction. Since q in (|27|) has compact 
support in Jl, the scattering transform depends only on the traces on d^l of the scattering solutions 
and ^i. Let = z^i + iz/2 be the complex-normal to the boundary. Then 

t(k) = -- [ e^q{z){'i'riz,k) -i^i{z,k))dn{z) = - [ e^ (d^{z, k) - id¥i{z, k)) dn{z) 
Jn ^ Jn 



2^ 



e''''v{z) {^r{z, k) - i^i{z, k)) da{z), (30) 
an 
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The last equality uses the fact that d{e'^'') = . 

Next we show how to reconstruct traces of and to dfl from the Dirichlet to Neumann 
map Afej^feg- First we reconstruct traces of Wr and Wi to (9$!. 

As in Knudsen and Tamasan we consider the single layer potential operator Sk ■ C°(f?0) — > 
a = l-2/p, defined by 

s'kfiz) = ^p.v. [ /(CMC - z)dC, z e dn, 

2vr« Jen 

where gk{z) = e~*^^/(7rz) is a Cauchy kernel for d which also takes into account the exponential 
behavior at infinity. For Lipschitz boundary is a bounded operator (e.g. see Muskhelishvili jl5j). 
Since q has compact support we have that Wr and Wi are analytic outside O and behaves like e^^^ 
at infinity. Traces of such functions will satisfy a singular boundary equations involving Sk- Inside 
they satisfy a pseudo-analytic equation. This will impose constrains (in terms of At^^fcj) their 
trace. We will prove that these two conditions are sufficient to determine the traces. 

We notices already that Wr{-, k),Wi{-, k) G C"(M^) with a = 1 — 2/p, whence their traces on 
dn are in C"((90). Let 



C^idn) := {h G C^idn) : / h{s)ds = 0}. 

JdQ 



Ian 

Define now a right inverse of the tangential vector field dg (here s is the arc length) on dfl by 

d;'f{t)= f f{s)ds, (31) 
Jo 

for / € C"(c?r2). In the above integral we fixed an arbitrary point on dO, from where we measure the 
arc length counter-clockwise. Notice that d^^ : Cg (50) — s- C^"^°(50) is a well defined (independent 
of the reference point) bounded operator. The following result defines a Hilbert transform for the 
pseudo-analytic maps. 

Lemma 3.4. Hb = —Ab-^^b2d'^^ '■ C'o (^^) ~^ C°'{dQ) is a bounded operator. 

Proof. Let g = d'^ f € C"+i(90) C W'^'^/P'P. Classical theory of PDE (e.g. see Gilbarg and 
Trudinger [S]) gives that the boundary value problem 

du du 

Au{x) +bi — (x) + 62 (x) — (x) = 0, xen (32) 
ox ox 

u\dnix) = g{x), X e 

has a unique solution up to a constant in Ty^'P(J7) and ||m||vk2.p(c) ^ C'||g||^2-i/j5,,5(gf^). Using the 
mapping properties of the Dirichlet to Neumann map we have 

\\T~(-bf\\c"{dn) < \\^bi,b29\\w^'P{dn) < l|V^lliyi-p(f7) ^ Il'"llw2>p(n) 

< C\\g\\w^-i/p,P(^QQ) < C\\g\\ci+'^(dn) < C'll/llc"(90)- 

□ 

Next we show that Tib reconstructs traces of the exponentially growing solutions on 90. 
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Theorem 3.5 (Trace theorem). Let b S Wc'^{i^). Consider the class of functions 

B = {he C^idfl) : lm{uh) e C^idfl)}. 

Then, for each A: G C arbitrarily fixed, the traces hr = Wr{-, k)\Q^., respectively hi = Wi{-, ^""6 
the unique solution in B of the systems 

(I - iSk)hr{z) = 2e''^, z G (33) 
Hi,{lxn{vhr)){z) = Re(i//ir)(z), z e dQ., (34) 

respectively, 

{I - iSk)hi{z) = 2ie'^^ , z e dfl, 
nb{lni{L'hi)){z) = Re{l^h^){z), z € dn. 

Proof. We argue only for Wr, the arguments for Wi are similar. 

We prove first the necessity. The arguments for ()33|) are identical to the ones in jl3j reason for 
which we only sketch them here. Fix A; G C and suppress the k dependence, we have Wr{-) = Wr{-, k) 
is analytic outside and e~'^^^Wr — 1 G LP{C — Cl). The Green-Gauss formula for z G C — gives 

,-..H^;m-i = --L/ £ri!l!yO,c^ (35) 

2vr2 Jen C- z 

Now let z approach (from the exterior) a boundary point zq and use Plemelj formula (see Muskhe- 
lishvih [13). 

hm — / = --f{zo) + —p.v. / 

to get (|33jl . 

Next we prove the necessity of H34() . Recall from lemma that Wr{z) = du{z) for some 
u G PF^'^(r2) which solve the equation (|24|) . Therefore 

hr = Wr\dQ. = ]j^{dx + idy)u\gQ. (36) 

For z G do, let {i'i{z),i'2{z)) be the unit outer normal, we also let i/(z) = I'l^z) + ii'2iz)- Next 
we express the partial derivatives for points on the boundary in terms of the tangent ds and the 
normal di, derivatives 

Vuix)=(-^' "^Va^^"' y (37) 

where we recall d^u = Ab-^^^h^u. Therefore 2h,. = {dx — idy)u = —iVdsU + VA.b^^i,,^u, or, using vV = 1, 

2vhr = Afo^^bjii — idsU. (38) 

Note that Im(i//ir) = —dsu/2 and thus hr (z B and (?7"'^(Im(i^/ir.)) makes perfect sense. Identifying 
the real part in gives H34() . Notice not only that we proved necessity but also we provided 
existence of solutions for (|3H) and (|33|) . 

Conversely, let /i G S be a solution of the system (|34|) and (|33)) . We extend h inside Vt by the 
following procedure. Inspired by (|38() define 5 = — Im(2z//i) G C'^{dQ.) then uniquely solve the 
boundary value problem © for u G VF^'^(r2). Notice 5 is real valued hence u has also real values. 
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Define Wr{z) = du{z) inside and notice that du\Qfi G (7"((?r2). Now check that OuIqq = h. 
Indeed, as before, 29n = —iVdsU + VKi,-^^h^u = iVlm{2vh) + I7Re(2z^/i). The last equahty used the 
fact that his a. solution of (|34() . Multiplication by v gives du = h. 
Inspired by (|35j) define Wr analytically outside Q by 



WJz) = e'''' / z G C - 0. (39) 

271-1 Jan C,- z 

The fact that h solves (|33() implies that lim^^^gggf^ Wr{z) = h{zQ). Thus Wr is an outside continuous 
extension of h. Moreover, e~'^^^Wr — 1 = 0{l/z) for z large, hence Wr G L^(C — il). 

We produced a continuous map in which solves H25|) both inside and outside $7 and behaves 
like e*^'^ for z large. We need to check that it solves the equation (|25)) across the boundary. Since b 
has compact support inside O we have that Wr is in fact analytic in both sides of the boundary and 
continuous across. Morera's theorem asserts that Wr must be in fact analytic across. Therefore 
Wr solves ()25() in the whole plane and has the right behavior at infinity. Uniqueness in lemma 
concludes the proof. □ 

Immediate consequence to the proposition 13.31 and to the pointwise estimates (|29j) we can 
determine the traces on dVt of ^r and ^i. Moreover by formula (|3U() we determine the scattering 
transform. 

Corollary 3.6 (Reconstruction of the scattering transform). Under the assumptions of the 
proposition we have 



e"^ ^{z) = \\mWr{z,k), z € dQ. (40) 

fe^oo 



and for any k C we recover 



^r{z, k) = e^ ^\z)Wr{z, k), ze on, (41) 
^i{z, k) = e^'^^{z)Wi{z, k), ze dn. (42) 



Moreover, 



m = ^ I e'^'^Viz) [^r{z, k) - i-^i{z, k)) daiz), (43) 
Jan 

is a function in L''(R^) n L^(M^) n L'''(M^) for some r < 2, r^^ = r^^ - 1/2 and r"^ + r'^^ = 1. 

Now we use the inverse scattering method of theorem 1.2 to reconstruct q. 
Corollary 3.7. Let ^r ~ e*^^ and ~ ie^^'' in U' for large k € C be the unique solutions 

^{k) + t{k)'^{k) =0, k£C. (44) 
dk 

Then 

e^t{k){^r{z,k) -i<^i{z,k))dfi{k). (45) 



Knowing q we also know |6| since from ()27() we have \q\ = \b\. Next we show how to determine 
its argument by solving (|27jl to recover b. 

The following result is due to Cheng and Yamamoto |5] . For the sake of completeness we sketch 
its proof. 
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Lemma 3.8. If q ^ LP{Q) then there exist at most one solution b E L^{^) of the equation 

b{z)e^''^^'^-^''^^'^ = q{z), zen. (46) 

Proof Assume there are two solutions 61, 62 G LP{Q.) and let d = d^^{b2 - 61) G W^'P{Q) C H'^{Q). 

From (|4U() we have d = on dQ. Hence d G Hq{Q). Since both solve ()46p we have bi{z) = 
b2{z)e'^^^^~^^^\ from where 

\dd\ = I62I • le'^^^)-^^") - 1| = \q\ ■ \e^(^)-di^) - l\ < \q\ ■\d-d\< 2\q\\d\. (47) 

By Carleman estimates for d G H^(M?) of compact support (see Hormander Prop. 17.2.3) 
we have 

[ A^\d\'^e^^dx < 4 / ladpe^^dx, (48) 

for some strictly convex function (/? G (7^(17). Approximate a (/? G H'^'^\Q) solution of A(/? = 17|gp 
in L^(i7) by a smooth sequence ip^ ^ uniformly on $7. Then 

/" Av9„[dpe2'^"dx < 16 /" |gp|(i|2e2'^"dx. 

For n sufficiently large the reverse inequality holds. Hence d = and 61 = 62. □ 

We are left to find the unique solution of 

Lemma 3.9 (Phase unwrapping). Let u G 1 + L^(M^) be the unique solution of 

dv = qv (49) 

then V vanishes on a set of measure zero. Define b = qv/v on the set where v does not vanish, else 
we can set b = q. Then b is the unique solution of ^ in LP{n). 

Proof. Existence and uniqueness of v follows from the Fredholm alternative as before. It is known 
from Vekua |V62j that the set of zeroes of pseudo-analytic functions has measure zero . Since 

bv = qv we have that v also solves dv = bv in the whole plane. Equivalently d{e~^ ^v) = 0. Thus 
g-9 by jg analytic and also goes to 1 as \z\ — > 00. By Liouville's theorem we have v = ^ . From 
its definition we have 

□ 

4 Concluding Remarks 

In order to solve the inverse problem one finds first the traces of the exponentially behaving solu- 
tions. It is easy to show that any solution of outside a disk can be represented as a series 

W{z.k)=e-^Y.Vn^ 

n=0 

with a„'s unknown coefficients. We determine them by solving the singular boundary integral 
equations (|34|) . This step is severely ill posed and regularization techniques are necessary, truncation 
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in the above series helps, see Knudsen ^1] for further ideas of regularization. Moreover, there is 
only a logarithmic type stability, see Barcelo et. al. in J]. 

The second step consists in constructing the scattering transform t{k) via the formulae of 
corollary ESI Next we solve the weakly singular integral equations Q in the fc-space. This part is 
stable. It is here that we need the e-extra regularity. One needs t S L'"(M|) for some r < 2 in order 
to solve If g is only in Lc then t £ L^(M^) (according to Sung [53] as corrected by Brown and 
Uhlmann 4 ) and this suffices for uniqueness. This would recover the uniqueness result of Cheng 
and Yamamoto. It is not clear how to find solutions of © when t G L^(M^). 

Reconstruct q from the formula 0. Notice that we have estimates of decay in k for t € L''(M|) 
as well as for e~^^^{^r — i^i) — 2 as given in Q). These can lead to estimates of the truncation 
error in the integral in (|7|). 

One of the questions in concerned the characterization of traces of exponentially behaving 
solutions in the first order system in Q: dv = qw and dw = qv. A partial answer was given by 
Knudsen and the author in ^3] for q of the special form q = df with / real valued. We can give 
now the answer for a general q. Note that v ±w solves the 5-equation du + qu = 0, and that we 
characterized the traces on the boundary of such solutions in terms of a Hilbert transform. 

Acknowledgement 

I would like to thank Professor A. Nachman for his generous sharing of ideas in the inverse scattering 
theory, and the organizers. Professors J. McLaughlin and H. Engl, for the invitation in the special 
semester on Inverse Problems at IPAM-UCLA. 

References 

[1] J. A. Barcelo, T. Barcelo and A. Ruiz 2001 Stability for the conductivity equation in the plane 
for less regular conductivities, J. Differential Equations, 173 231 - 270 

[2] R. Beals and R. R. Coifman, 1988 The spectral problem for the Davey-Stewartson and Ishi- 
mori hierarchies. Nonlinear evolution equations: Integrability and spectral methods, Manchester 
University Press, Manchester, 15-23 

[3] R. M. Brown and G. Uhlmann, 1996 Uniqueness in the inverse conductivity problem with less 
regular conductivities in two dimensions. Private notes 

[4] R. M. Brown and G. Uhlmann, 1997 Uniqueness in the inverse conductivity problem for 
nonsmooth conductivities in two dimensions. Comm. Partial Differential Equations 22 1009- 
1027 

[5] J. Cheng and M. Yamamoto 1998 Unique Determination of two convection coefficients from 
Dirichlet to Neumann map. University of Tokyo preprint UTMS 98-31, to appear in SIAM J. 
Math. Anal. 

[6] J. Cheng and M. Yamamoto 2000 The global uniqueness for determining two convection coef- 
ficients from Dirichlet to Neumann map in two dimensions Inverse Problems 16 L25-L35 

[7] C.Tong J. Cheng and M. Yamamoto 2004 Reconstruction of convection coefficients of an elliptic 
equation in the plane by the Dirichlet to Neumann map, electronic preprint arxiv.org 

[8] R. R. Coifman and Y. Meyer 1978 An deld des operateurs pseudodifferentiels, Asterisque 57, 
Societe Mathematique de Prance, Paris 



12 



[9] D. Gilbarg and N. Trudinger 1983 Elliptic partial differential equations of second order, 
Springer- Verlag, Berlin 

[10] G. M. Henkin and R. G. Novikov 1987 5-equation in the multidimensional inverse scattering 
problem, Russian Math. Surveys 42 109-180 

[11] L. Hormander 1985 The Analysis of Linear Partial Differential Operators, vol 3, Springer 

[12] K. Knudsen and A. Tamasan 2001 Reconstruction of less regular conductivities in the plane, 
MSRI preprint 2001-035. 

[13] K. Knudsen and A. Tamasan 2003 Reconstruction of less regular conductivities in the plane, 
MaPhySto preprint 1398-2699, in print in Comm. Partial Differential Equations 29 

[14] K. Knudsen 2003 A new direct method for reconstructing isotropic conductivities in the plane, 
Physiol. Meas. 24 391-401 

[15] M. I. Muskhelishvili 1953 Singular integral equations. Boundary problems of function theory 
and their application to mathematical physics, P. Noordhoff N. V., Groningen 

[16] A. I. Nachman, M. J. Ablowitz 1984 A multidimensional inverse-scattering method, Stud. 
Appl. Math. 71 243-250 and 251-262. 

[17] A. I. Nachman 1988 Reconstructions from boundary measurements Ann. of Math. 128 531-576 

[18] A. I. Nachman 1996 Global uniqueness for a two-dimensional inverse boundary value problem, 
Ann. of Math. 143 71-96 

[19] L. Nirenberg and H. F. Walker 1973 The null spaces of elliptic partial differential operators in 
M" J. Math. Anal. Appl. 42 271-301. 

[20] E. M. Stein 1970 Singular integrals and differentiability properties of functions, Princeton 
University Press, Princeton, N.J. 

[21] J. Sylvester and G. Uhlmann 1987 A global uniqueness theorem for an inverse boundary value 
problem, Ann. of Math. 125 153-169 

[22] L. Sung 1994 An inverse scattering transform for the Davey-Stewartson II equations. I J. Math. 
Anal. Appl. 183 121-154 

[23] L. Sung 1994 An inverse scattering transform for the Davey-Stewartson II equations. II, J. 
Math. Anal. Appl. 183 289-325 

[24] L. Sung 1994 An inverse scattering transform for the Davey-Stewartson II equations. Ill J. 
Math. Anal. Appl. 183 477-494 

[25] G. Uhlmann 2003 Inverse Boundary Problems in two Dimensions Function Spaces, Differential 
Operators and Nonlinear Analysis- The Hans Triebel Anniversary Volume, eds. D. Haroske et 
al. Birkhauser, Basel-Boston-Berlin 183 - 203 

[V62] I.N. Vekua 1962 Generalized analytic functions, Pergamon Press, London 



13 



